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The concertina is a magnetization pattern in elongated thin-film elements of a soft material. It 
is a ubiquitous domain pattern that occurs in the process of magnetization reversal in direction of 
the long axis of the small element. 

Van den Berg argued that this pattern grows out of the flux closure domains as the external field 
is reduced. Based on experimental observations and theory, we argue that in sufficiently elongated 
thin- film elements, the concertina pattern rather bifurcates from an oscillatory buckling mode. 

Using a reduced model derived by asymptotic analysis and investigated by numerical simulation, 
we quantitatively predict the average period of the concertina pattern and qualitatively predict its 
hysteresis. In particular, we argue that the experimentally observed coarsening of the concertina 
pattern is due to secondary bifurcations related to an Eckhaus instability. 

We also link the concertina pattern to the magnetization ripple and discuss the effect of a weak 
(crystalline or induced) anisotropy. 



I. INTRODUCTION 

To our knowledge, the te rm conce rtina was introduced 
by van den Berg et. al. in |vdBV82l p. 880]. In that pa- 
per he explains the formation of this domain pattern in 
Permalloy thin-film elements that are fairly thick (thick- 
ness t = 350nm), with a rectangular cross section (width 
i = 15|LLm) that is not too elongated (length 60|LLm). After 
near-saturation along the long axis, a concertina pattern 
grows out of the flux closure domains at the short edges of 
the cross section during subsequent reduction of the ex- 
ternal field i^ext 7 until the pattern eventually invades the 
entire sample, see Figure [T] on the right. Our experimen- 
tal observations suggest that in very elongated samples 
(length 2nm, thickness 10 to 150nm, width 10 to 100|LLm) 
a bifurcation is at the origin of the concertina. As a con- 
sequence, the pattern forms simultaneously all over the 
sample. We will report on van den Berg's explanation in 
more detail after introducing the micromagnetic model. 



A. The micromagnetic energy 

The variation of the applied magnetic field in the ex- 
periments is on a very slow time-scale so that the magne- 
tization always relaxes to equilibrium. Hence our theoret- 
ical analysis is based on the micromagnetic (free) energy, 
which we introduce now. The magnetization of a ferro- 




FIG. 1: Concertina in a very elongated (length 2 mm) sample 
of width 50 i^m and thickness 50 nm (left) and in a sample of 
width 35 p.m, thickness 40 nm and moderate length 110 ^im 
(right). The left image shows only the center of the stripe 
which is less than 10 percent of the whole sample. As indi- 
cated by the blue arrows, the gray scales encode the transver- 
sal component of the magnetization in the domains. 



magnetic sample occupying some domain Q is described 
by a vector field m = (mi, 7712,7713). The micromagnetic 
energy £^(777,) is given by 

E{m) =d^ f |Vm|2 dx+ / li^^strayl' dx 

J Q, J all space 

- Q / {m-ef dx-2 [ 



(1) 



dx — 2 I i^ext • ^ dx. 



The energy ([T]) is already partially, i.e., except for 
lengths, non-dimensionalized. In particular the magneti- 
zation has length one. Outside of the sample, it vanishes 
identically: 



1 in the sample 



m = outside of the sample Cl. 



(2) 



2 




FIG. 2: Formation of the concertina pattern in the experiment: The pictures show a section near the center of the elongated 
thin- film element. The unstable mode grows into a domain- wall pattern which coarsens several times. The two upper series 
show a sample of 30nm thickness of low anisotropy. The two lower series show a sample of 30nm thickness of higher anisotropy. 
The width is 30|J,m and 50|J,m, respectively. 



Let us now briefly introduce and discuss the different 
energy contributions: 

The first contribution in ([T]) is the so called exchange 
energy. (The gradient acts component wise, i.e., |Vmp = 
X^i^i This term favors a uniform mag- 

netization. The material parameter d is called the ex- 
change length and measures the relative strength of ex- 
change with respect to stray-field energy, see below. The 
exchange length is typically of the order of a few nano- 
meters. 

The second contribution in ([T]) is the stray- field energy. 
The static Maxwell equations state that the magnetiza- 
tion m generates a stray field i^stray that is described 

by 



V X i^stray — 

V • (i^stray + m) = 0, 



(3) 



where both equations hold in the whole space and B — 
^stray + m is the magnetic induction. Hence the stray 
field is generated by the divergence of the magnetization. 
Since the magnetization is discontinuous at the boundary 
of the sample dVt^ cf. ([2]), the second equation in (j3]) has 
to be understood in the following sense: 



^ ' -f^stray 



— V • m in the sample 
outside of the sample 1], (4) 

[^stray ' y\ = vfi • V ovi the bouudary 91], 

where v denotes the normal to the boundary of the sam- 
ple, and [i^stray • ^] denotes the jump that i^stray • ^ ex- 
periences across the boundary dVl. Hence we distinguish 
two different types of sources of the stray field - in anal- 
ogy to electrostatics one commonly speaks of charges - 
namely 

magnetic volume charges V • m in and 
magnetic surface charges v - m on dVt. 

Later on we will also use the following equivalent distri- 
butional formulation of ([4]), which is obtained by testing 
with smooth functions ^ vanishing at infinity, namely 



all space 



stray 



VC dx 



The third contribution in ([T|) models a uniaxial 
anisotropy, i.e., the preference of an easy axis e = 
(61,62,63) in a material. The material parameter Q > 
is called the quality factor. It measures the relative 
strength of anisotropy with respect to stray-field energy. 
A uniaxial anisotropy can for example come in form of 
crystalline or induced anisotropy. Notice that the poly- 
crystalline anisotropy in a material like Permalloy can be 
described as a position-dependent easy axis e(x). 

The last contribution in ([T]) is called Zeeman energy. 
This term models the interaction and favors the align- 
ment of the magnetization with an external magnetic 
field i^ext- 

The specific material parameters of our samples are 
discussed in Section [TdI 



B. Van den Berg's explanation of the concertina 



m • VC dx. (5) 



We now return to van den Berg's paper |vdBV82j . 
Combining the explanation of the formation of t he con- 
certina p attern therein with the insight from |BS89L 
we can give the following updated version of 
the explanation presented in [vdBV82, Sections A & B]: 
It is a theory on the two-dimensional mesoscopic mag- 
netization pattern; by two-dimensional, we understand 
that the magnetization is in-plane, i.e., ms = 0, and inde- 
pendent of the thickness direction, i.e., m = m(xi,X2); 
by mesoscopic, we understand that the walls are replaced 
by sharp discontinuity curves that are charge-free in the 
sense that the normal component of the magnetization 
does not jump; moreover, the magnetization is tangential 
to the lateral edges of the sample so that there are no sur- 
face charges. In sufficiently large thin- film elements and 
for sufficiently low external fields, [BS89J now postulate 
that the two-dimensional mesoscopic magnetization pat- 
tern arranges itself in such a way that the corresponding 
continuous magnetic charge density a = —{dimi -\-d2m2) 
generates a stray field Hgtray that expels the external 
field i^ext from inside of the sample (like in electrostat- 
ics). 

In jDKMO05l |. see |DKM+Olj for an efficient account, 
it is shown that in the regime of sufficiently large thin- film 
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elements (i.e., t <^ i and it log | with comparable 

lateral dimensions of the order ~ this principle ex- 
tends to moderately large fields (of the order ~ |): In this 
case, the stray field i^stray in general can no longer expel 
the external field i^ext everywhere in the sample, since 
the (total) charge density a = — {dim 821712) is limited 
hy mf -\- ml = 1. The charge density a is uniquely de- 
termined by a convex variational problem only involving 
the stray-field energy and the Zeeman energy. At least 
some aspects of the mesoscopic two-dimensional magne- 
tization pattern (mi, 7712) can be recovered from a: The 
characteristics of (mi, 777-2), i-e., the curves a long whic h 
(mi,m2) is normal (called "trajectories in |vdBV82| ). 
have curvature given by a. However, due to the poten- 
tial discontinuity curves of the mesoscopic magnetization 
(mi,m2), this seemingly rigid condition does not suffice 
to determine (mi,m2) - even the fact that the discon- 
tinuity curves are charge-free is still not enough. Notice 
that it is easy to construct a particular solution (mi, m2) 
for any charge density a v ia the max imal solution of a 
modified eikonal equation jPKM+QlL p.2987]). On the 
other hand, in the region where the external field has 
penetrated, t he magnet ization (mi,m2) is uniquely de- 
termined , cf. IDKM+Qll . p.2987], and has no discontinu- 
ities, cf. jvdBV82L D.883]. 

Van den Berg gives a recipe how to construct a so- 
lution that corresponds to the experimental observation 
of a concertina pattern growing out of the fiux closure 
pattern: For sufficiently large external fields (i^ext = 

(/7ext,0,0), /lext > |), ^ext H stray doCS UOt vauish 

in the sample besides in the vicinity of the two distant 
edges; as a consequence walls only occur in the two fiux 
closure pattern there. As the external field is reduced, the 
penetrated region shrinks as the walls invade the sample. 
Each of the two fiux closure pattern has a "doublet" , that 
is, a point on one of the long edges where two wall seg- 
ments intersect. Note that the doublets were created as 
the central 180°-wall of the initial Landau state moved 
towards the edge where it broke up due to the apphca- 
tion of a large external field at the very beginning. The 
inner, i.e., most distant to the short edges, ones of the 
doublet walls fade out in the middle (with respect to the 
long edges) of the cross section. Van den Berg postulates 
that the position of the doublets does not change as the 
external field decreases further. As a consequence, each 
of the two inner walls grows - necessarily in direction of 
the characteristic - till is hits the opposite edge. There 
it must generate a "triplet" (a point on the edge where 
three walls meet); the middle wall must coincide with the 
previous one originating in the doublet. Again, as the ex- 
ternal field is further reduced, the position of each triplet 
is supposed to be fixed, the inner of the three walls grows 
towards the opposite edge. 

This process repeats itself till the two half-concertina 
structures growing from the short edges are linked in the 
middle (with respect to the short edges) of the cross sec- 
tion. For very elongated samples of length L ^ £^ the 
linking is expected at a field strength of order /lext ^ 




FIG. 3: The creation of a triplet out of the initial do ublet as 
described and sketched by van den Berg in jvdBV82j 

t£L~^ \nt£~^ and thus differs from the field at the be- 
ginning of the growth process by a factor (up 
to a logarithm). Speaking in mathematical terms, van 
den Berg postulates that the positions of doublets and 
triplets remain fixed as the field is decreased and appeals 
to continuity, i.e., the pattern should depend continu- 
ously on the value of the external field, to overcome the 
non-uniqueness of (mi,m2) mentioned in the previous 
paragraph. 

C. Van den Berg's vs. our explanation 

Our explanation for the formation of the concertina 
pattern is very different from the one of van den Berg. 
Instead of a successive outgrowth (along the sample) of 
the closure domains, we explain the concertina as a si- 
multaneous outgrowth (along the sample) of an unstable 
mode, best captured in very elongated thin-film elements. 
Indeed, our experiments were performed on thin-film el- 
ements of thicknesses t in the range of lOnm to 150nm, 
widths i in the range of 10|LLm to 100|LLm, but lengths in 
the range of 2mm. We recorded the pattern at three dif- 
ferent sections, equidistant and equidistant to the short 
edges of the cross section, and observed qualitatively the 
same pattern at the same values of the external field. 

Not surprisingly, our theoretical predictions are quite 
different from those in |vdBV82| - aheady in terms of 
scaling. Van den Berg's explanation entails two different 
scales of the external field 

• ^ex?" ^ I beginning of the build-up process 
and 

• hl^^ ~ t£L~'^lnt£~^ for the completion when the 
external field is totally expelled from the sample 

whereas in our case there is one critical field h^^^ at which 
the simultaneous formation of the concertina along the 
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sample - independently of the specific position - due to 
an interior instability begins. This critical field is given 
by K^t ^ -d^l'^t'^l'^t^l'^, see Regime III in Subsection 
II F[ for isotropic samples - thus the instability would 
only occur after the field is reversed and thus when the 
van den Berg concertina has already invaded the sam- 
ple. However, as we discuss later, the critical field is 
shifted in case of a transversal anisotropy /i*^^ ^ext+Q- 
It turns out that even for relatively weak transversal 
anisotropy the formation thus starts before the field is 
reversed, see Section [VH a). For very elongated samples, 
i.e., L > ^, we have that /i^^f^"" > and = 
for the limit case of an infinitely extended sample. The 
strength of the anisotropy and the geometry of the ma- 
jority of the samples that we investigated is such that 
^exf > ^ext > ^ext^- We thus cxpcct the following sce- 
nario in very elongated samples: At the van den 
Berg build-up process starts at the tips of the sample. 
As the field is reduced, the concertina grows slowly into 
the sample from the tips. Meanwhile, as h%^^ is attained 
our instability occurs all over the sample - sufficiently 
far away from the tips and way before the van den Berg 
Hnking could take place in the center of the sample, see 
Table H 





7 begin 
'''ext 


^ext 


7 end 
'''ext 


Weak anisotropy 
(Q = 1.3 X 10"^) 








i = 30|^m 


1 X 10"'' 


-4.0 X 10""^ 


1.6 X lO"*" 


I = 50M,m 


6 X lO"'' 


4.2 X 10"^ 


2.8 X 10""* 


Stronger anisotropy 
(Q = 5.0 X 10"^) 










1 X 10-^ 


3.4 X 10-"^ 


1.6 X 10-^ 


i 50M,m 


6 X 10"^ 


4.2 X lO"'' 


2.8 X 10"^ 



TABLE I: Comparison of the characteristic fields in van den 
Berg's theory of the concertina and the critical field in our 
instability for the samples shown in Figure [2] (thickness t — 
30nm, length L — 2mm. Apart from the sample of weak 
anisotropy and small width, the characteristic fields appear 
in the expected order. 

Whereas in |vdBV82| the appropriate scale for the con- 
certina width w is given by ^ - in particular independent 
on the thickness t - it is given by and dP'^'^P^'^t~^^'^ in 
our case, in qualitative accordance with our experimental 
observations illustrated in Figure [H 




FIG. 4: Concertina in Permalloy samples of width £ — lOO^im 
and thickness t — 30nm (left), t — 80nm (center), and t — 
300nm (right). The width average period of the pattern is a 
decreasing function of the thickness t. 



D. Experimental setup and samples 

In the experiments, we investigated magnetic films 
of nano-crystalline Permalloy (NigiFeig) and amorphous 
Co6oFe2oB2o of various thicknesses and induced magnetic 
anisotropy values. The films were deposited by mag- 
netron sputtering under ultra-high vacuum conditions. 
In order to control the grain growth of the polycrys- 
talline films, a Ta seed (5 nm) layer was used for the 
NigiFeig deposition. In all cases a magnetic in-plane sat- 
uration field was applied during film deposition to control 
strength and direction of the induced anisotropy. Using 
different magnetic field histories, films of different effec- 
tive induced anisotropy were obtained. 

• A first set of samples was deposited in the presence 
of a homogeneous, static magnetic field. This re- 
sults in a maximal and also well-aligned induced 
uniaxial anisotropy. A series of Permalloy and 
CoFeB samples was obtained by this method. 

• In a second set of samples the induced anisotropy 
was strongly reduced. In order to ensure this, the 
films were deposited in a magnetic field of alter- 
nating orthogonal alignment. The field direction 
was changed after approximately every 5nm of film 
growth. The superposition of the so-obtained or- 
thogonal anisotropy axes results in a strongly re- 
duced induced anisotropy. 

The relevant material parameters - for the comparison 
of the experimental observations to the theoretical pre- 
dictions - are the following: 

• Exchange length d: Permalloy 5nm, CoFeB 3nm. 

• For both materials the saturation polarization is 

~ IT and the stray- field energy density is given 
hy Kd-4.x 10^ J/m^. 

• The uniaxial anisotropy coefficient is j^Permaiioy ^ 
200 J/m^ for the high- anisotropy Permalloy and 
^CoFeB ^ gQQ CoFeB, respectively. 
For the low- anisotropy Permalloy films we have 

^Permalloy ^ 5Q 

• Quality factor Q = Ku/Kd'- High and low- 
anisotropy Permalloy 0.5 x 10~^ and 0.125 x 10~^, 
respectively, and CoFeB 1.5 x 10~^. 

• The average size of the individual grains of Permal- 
loy is ^grain ^ 12 to 15nm. It is assumed that up to 
a film thickness of about 30nm, the grains display 
a column-like shape. 

• The film thicknesses range from 10 to 150nm, the 
film widths from 10 to 100|LLm. 

After film deposition, elongated stripes of various widths 
and a length of 2000 |LLm were patterned by photolithog- 
raphy and subsequent ion beam etching. The stripes 



5 



were aligned, both, parallel and orthogonal to the in- 
duced anisotropy direction. 

The observation of domains and magnetization pro- 
cesses was carried out in a digitally enhanced Kerr micro- 
scope, ^HS98:] . The longitudinal Kerr effect was applied 
with its magneto-optical sensitivity direction transversal 
to the stripe axis. The dominant wavelength of the ob- 
served concertina patterns was computed by Fast Fourier 
Transform. The result of the computation is in agree- 
ment with the average wavelength determined by man- 
ually counting the folds in the images, as soon as the 
concertina becomes discernible to the eye during field 
reduction. The typical strength of the magnetic field, 
which is applied for the saturation at the very beginning, 
is of the order of some mT. 



E. Nucleation 

We are interested in the magnetization pattern in elon- 
gated thin-film elements of width £ (in X2 -direction) and 
thickness t <C ^ (in xs-direction) that forms under the 
variation of an external field aligned with the long axis 
(the xi-axis), that is, of the form i^ext = (— ^ext,0,0). 
(To simplify the notations, the minus is introduced which 
entails a positive critical field, see below.) As indi- 
cated above, we observe no infiuence of the sample's 
short edges on the formation of the concertina away from 
the short edges. Since it greatly simplifies the theo- 
retical treatment, we therefore henceforth will assume 
that the sample is infinite in xi-direction (and occasion- 
ally, for instance in the numerical treatment, impose a 
large, but artificial period in that direction). One conse- 
quence of that assumption is that the uniform magneti- 
zation m* = (1,0,0) is a stationary point of the energy 
functional - that is, satisfies the corresponding Euler- 
Lagrange equations that express a torque balance at ev- 
ery point of the sample - for all values /lext of the external 
field i^ext of the form above. 



X2 




Xi 



FIG. 5: The idealized geometry of the sample. The homo- 
geneous external field ifext is parallel to the long axis of the 
sample. 

The experiments suggest that as the strength of the 
field is reduced starting from saturation, i.e., he^t < 0, 
and finally reversed, a bifurcation at some critical value 
^ext > of the external field i^ext = (— ^ext, 0, 0) is at the 
origin of the concertina pattern, see Figure [2] and Figure 
[4Ql Due to the unit-length constraint ([2|), infinitesimal 
variations of m* are of the form Sm = {0,6m2,Sms). 
Since the uniform magnetization only generates Zeeman 
energy, the linearization of the energy in m* - neglecting 
anisotropy - is given by the exchange energy and the 



stray-field energy of the infinitesimal variation, that is. 
In I^^^P + /all space l^stray ((^^) P dx, augmented 
by the linearization of the Zeeman energy. The latter is 
due to the constraint ([2]) given by — /lext /^(^^2 +^^3)x 
which is a consequence of the expansion mi = (1 — Sml — 

F. Unstable modes 

We start with the linear stability analysis of the uni- 
form magnetization by discussing potentially unstable 
modes on the level of the linearization of the energy. We 
argue at which value of the external field /lext each of the 
modes becomes unstable. At this so-called critical field 
the infinitesimal release of Zeeman energy becomes larger 
than the infinitesimal contributions due to exchange and 
stray-field energy. We neglect uniaxial anisotropy (i.e., 
we set Q = 0) for the moment, since on the level of 
this infinitesimal discussion, a longitudinal or transver- 
sal anisotropy just leads to a shift of the critical field 
^ext ^ ^ext + Section [Vll Since the shift entails 

that the sign of the critical field can change, we note in 
that context that if we speak about reducing the strength 
of the external field we usually mean that the critical field 
is approached from saturation (/lext = — oo) if not stated 
differently. Similarly we say that the external field is in- 
creased after the critical field is passed. In this sense, the 
critical field is interpreted as the zero point on the scale 
of the external field, cf. Figure [6l 

m* stable m* unstable 



saturation /lext = — oo /igxt > ^ext 
FIG. 6: The scale of the external field /lext- 

I) The first mode we discuss is a coherent rotation, i.e., 
Sm = (0, (5m2, Sms) is constant in space. Such a mode re- 
leases Zeeman energy per length in xi -direction of the in- 
finitesimal amount he^t^tSA^^ where 6A = {Sm2-\-Sml)^^'^ 
is the infinitesimal amplitude of the coherent rotation. A 
coherent rotation necessarily generates surface charges. 
Since the top and bottom surfaces have larger area than 
the two lateral surfaces, an in-plane rotation {6ms = 0) 
is favored, cf. Figure [71 This mode generates surface 
charges of infinitesimal density ±M. Over distances i 
much larger than t, these surface charges act like two op- 
positely charged wires at distance i of line density tSA - 
also in the following if not mentioned otherwise always 
infinitesimally and per length in xi-direction. Hence 
the mode generates an infinitesimal stray field of order 
~ t"^ (in it~^)6A^ . Therefore, this mode becomes unstable 
when /lext ^ U-^{ln£t-^). 

II) The second mode we consider is buckling, cf. Figure 
[HI The magnetization avoids the lateral surface charges 
by just laterally buckling in the middle of the cross 



FIG. 7: Coherent rotation and generated surface charges. 



section, i.e., 6m = (0, Msin(7r^), 0). However, since 
V • Sm = 7r^~^Mcos(7r^), the surface charges of the 
coherent rotation turn into volume charges. At distances 
much larger than t from the cross section, these volume 
charges act like surface charges of amplitude ^ £~^t6A. 
Since these surface charges change sign over a distance 
they generate a stray field which extends a distance 
^ £ away from the cross section, and which is of the 
magnitude ~ i~^t6A. Hence this mode generates a stray 
field energy ~ t^3A^ ^ which is only smaller by a logarithm 
than in case of the previous mode of coherent rotation. 
Moreover, since \V6m\^ = tt^^^^M^ cos^(7r^), the mode 
generates exchange energy ~ d?£~^t8A^. Since the re- 
lease of Zeeman energy scales as ~ hQ^t^tSA^ as in case 
of the first mode above, this mode becomes unstable at 
^ext ^ in the regime t <C d'^i~-^ and at /lext ^ t£d~'^ 

in the regime t ^ 




FIG. 8: Buckling mode and generated volume charges and 
stray field - for reasons of a clear presentation only drawn in 
the region above the sample. 



III) The third mode we discuss is oscillatory buck- 
hng, cf. Figure [9l This mode reduces the stray-field 
energy through a modulation of the lateral buckling in 
xi-direction, i.e., Sm = (0, Msin(7r^) sin(27r^), 0) with 
a wavelength w that satisfies t <^ w <^ i. Since w ^ 
the volume charges generated by this mode act like sur- 
face charges of amplitude ^ £~^tSA over distances much 
larger than t from the cross section. However, these sur- 
face charges change sign over a distance <C ^, so that 
the generated stray field only extends over a distance ~ w 
away from the cross section. Hence this mode generates 
a stray- field energy ~ £~^t'^w5A'^, which is substantially 
less than the stray-field energy of the two prior modes 
for It; <C ^. However, since w <^ the exchange energy is 
now dominated by the oscillation in xi -direction, which 
leads to an infinitesimal exchange energy ~ (Piw~'^t6A^ . 
Hence the wavelength w which leads to the minimal in- 
finitesimal total stray- field energy and exchange energy of 
order - £/^£-^/H^/^SA^ is given by w* - d^/^^Vs^-i/s^ 
This is consistent with our assumption t <^ w <^ £ 
provided d'^£-^ <C t <C (d£Yl^. The oscillatory buck- 
ling mode becomes unstable at a field strength of order 

/lext ~ ^2/3^-4/3^2/3^ 

IV) The fourth mode we consider is curl- 



FIG. 9: Oscillatory buckling mode and generated surface 
charges and stray field - for reasons of a clear presentation 
only drawn in the region above the sample. 



ing. This mode avoids charges altogether by 
an X3-dependent magnetization, i.e., bm = 
(0 , M sin (tt ^ ) cos (tt ^ ) , 1 1 cos (tt ^ ) sin (tt ^ ) ) . The 
exchange energy in this case is dominated by the gradi- 
ent in X3-direction which scales as ~ d^£t~^3A^. Hence 
the curling mode becomes unstable at /lext ^ d^t~'^ . 

The discussion above shows that there are (at least) 
four different parameter regimes for the nucleation - ex- 
pressed in terms of the two non-dimensional parameters 
t/d <^ £/t. These regimes are characterized by a certain 
scaling of the critical field h^^^ in the sense that one of the 
modes becomes unstable as the external field passes the 
corresponding field, while the other three modes are still 
stable, cf. Figure [TOl In particular, the oscillatory buck- 
ling mode is the first mode to become unstable at a field 
hLt - d^l'^£-^lH^I^ in the regime d^£-^ <C t <C {d£fl^. 
By a refinement of the above discussion, it can be rig- 
orously shown that the re are exactly four regimes, cf. 
Theorem 1 in jCAO06al . p.357]. 



^\osc. buckling- IV 
- \ ■ curling 




FIG. 10: Phase diagram with the four regimes for the nucle- 
ation. 



G. Period of the unstable mode: Experiment vs. 
theory 

Clearly, the regime of interest for us is the Regime 
III, i.e., the oscillatory buckling regime characterized by 
d^£~^ <C t <C {d£Y^'^. In this regime, an asymptotic 
analysis of the linearization of the energy on the ba- 
sis of the above discussion shows that the (first) un- 
stable mode is indeed asymptotically of the f orm dm = 
(0,Msin(7r^)sin(27r^),0), cf. Theorem 1 in jCAOQGbl . 
p. 389] and see also below. Based on a refinement of the 
prior linear stability analysis one can moreover determine 
the asymptotic behavior of w* including the numerical 
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factor that is given by 

W* « (327r)l/3d2/3^2/3^-l/3_ (g) 

So far we have learned that in Regime III at field 
strengths h^^^ - d^l'^ i'^l'^t^l'^ there is a bifur- 
cation in direction of the unstable mode (5m = 
(0,Msin(7r^)sin(27r^),0). We claim that the con- 
certina pattern grows out of this unstable mode. If so, 
the experimentally observed period i^^exp should be close 
to the period it;* of the unstable mode. Defining and de- 
termining i^exp is delicate: As /lext increases (after the 
critical field /i*^^ is passed), there is a continuous tran- 
sition from the magnetization ripple, see Subsection [Vl 
to the concertina pattern, which is far from exactly pe- 
riodic, and which coarsens subsequently, see Section [TTTl 
As u'exp we take the average period as soon as the con- 
certina pattern is discernible to the eye. Figure [TTI shows 
the result of this comparison for a broad range of sample 
dimensions i and t and (therefore) a fairly broad range of 
periods it;*: The ratio of the smallest width t compared 
to the largest is 5, the ration of the smallest thickness t 
compared to the largest is 15. The smallest period it;* 
is expected for a thick film of small width, the largest 
period for a thin film of large width, differing by a fac- 
tor close to six (neglecting the prediction for the broken 

or defect samples). The ratio -^r- of the experimental 
period with respect to the prediction ranges around two. 
We basically see this as a confirmation of our hypothesis, 
namely that the concertina grows out of the oscillatory 
buckling. Notice that the deviation has a clear trend: 
i(;*xp is larger than it;*. We give an explanation for this 
systematic deviation in Section Hill 

H. A reduced energy functional 

In the forthcoming section we start with the investiga- 
tion of the type of the bifurcation. For the moment we 
continue to neglect anisotropy, although it may affect the 
type of bifurcation as we shall discuss in Section IVTl In 
order to understand the type of bifurcation we now first 
pass to a reduced model adapted to our Regime III: The 
form of the unstable mode suggests that the out-of-plane 
component and the dependence on the thickness variable 
are negligible, i.e., we assume ma = and m = m(xi, X2) 
respectively. Since the unstable mode varies faster in x\- 
direction than in X2-direction, we neglect |92mp with re- 
spect to l^imp in the exchange energy density. Since the 
oscillation in the sign of the charge density is on smaller 
length scales in xi-direction than in X2-direction, we ne- 
glect with respect to ln\ + ln\ in the stray-field energy 
density, where i^stray = (^1,^2,^3)- Finally, since we 
are interested in small deviations from m* = (1,0,0), 

we expand m\ = — vr3^ ~ 1 — so that we may 
neglect |Vmip with respect to |Vm2p in the exchange 
energy density. We also use this approximation in the 
charge density and in the Zeeman energy. Hence (up to 
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FIG. 11: The theoretical period of the unstable mode is in 
good correspondence to the measurements: The upper image 
shows the ratio of the experimentally observed period and the 
period of the unstable mode. The white patches correspond 
to broken or defect-ridden samples. The lower displays the 
ratio of the period w* and the smallest period, i.e., i6;*(^ = 
50^im,t = 150nm), at all. Both images share the same color 
map. 

an additive constant) we are left with the reduced energy 
functional 

E{rri2) ^ d'^t / {dirri2)'^ dxi dx2 
+ / (^1+^3) dxi dx2 dx3— /lext ^ / ^2 dxi dx2, 

J all space J Q' 

(7) 

where the stray field i^stray = (^1,0, /13) is determined 
via 

dshi -dihs = 0, 
/ {hidiC + hsdsC) dxi dx2 dxs 

J all space 

= t I {-^diC + m2d2C) dxi dx2 for ah C, 

(8) 

which is a consequence of the alternative formulation ([5]). 
Here, Cl' denotes the in-plane cross section of our sample 
n = n' X (0,t). 
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We note that the stray-field energy is only finite if m2 
vanishes at the lateral long edges, i.e., m2(xi, ^2) = for 
^2 = 0,^ (as is true for the unstable mode). Notice that 
(j8]) can be written as 

dihi + dshs = for xs ^ 0, 

2 

[^3] = t{-di^ + ^2^2) for Xs = 0, 

where [hs] denotes the jump hs experiences across xs = 0. 
This formulation shows that X2 is just a parameter in the 
equations for the stray field, which behaves like a two- 
dimensional stray field (in the xi^a-plane) generated by 

2 

the "line charge" t{—di^ -\- 821712). 

We note that the only non-quadratic term in the energy 

2 

comes from the non-linear charge distribution t{—di^-\- 
d2m2). This allows us to derive the scaling of the ampli- 
tude of the magnetization: It should be such that both 
terms in the charge distribution balance. Since in view of 
the unstable mode the typical xi -scale of the variations 
of 1712 is given by w"^ ~ d^/^-^^/^t"^/^, whereas the typical 

X2 -scale of variations of m2 is given by the sample width 

2 

^, the contributions di ^ and 92m2 balance provided the 
amplitude of 1712 scales as (P^^i~^^^t~^^^ . This suggests 
the following non-dimensionalization of length, and re- 
duced units for the stray field and the magnetization: 

m2 = ^2/3^-1/3^-1/3^^^ 

(9) 

After the reduction leading to ([7j) , only the stray- field de- 
pends on the X3-component for which the relevant length 
scale is of course the wavelength of the oscillation, cf. Sec- 
tion II F[ III and Figure [9] - the relevant length scale in 
case of the magnetization leading to the reduction ([7j) is 
of course the film thickness t. If we also rescale the ex- 
ternal field - in the same way as the stray field - and the 
energy itself according to 



we obtain the reduced rescaled energy functional 



(10) 

(11) 



E(m2) = I {811112) dxidx2 
/ (/ii + /13) dxi dx2 dx3 - /lext / rh\(lxidx2, 

J all space J Q.' 

(12) 

where the reduced rescaled stray field is determined by 
dihi + dshs = for xs ^ 0, 
[^3] = (-9i^ + 021712) for Xs 0, 



The reduced rescaled formulation shows that the reduced 
energy functional contains just one non-dimensional pa- 
rameter, namely the reduced external field /lext - in- 
stead of four parameters (exchange length, sample di- 
mensions and /lext) for the full model. Moreover, the 
vector ^e\d m = (mi, 7712,7713), function of three variables 
(xi, X2, X3), has been replaced by the scalar function 7712, 
function of two variables (xi,X2). Finally, the compu- 
tation of the stray field is a ^it'o-dimensional computa- 
tion (in {xi^xs) only with X2 as a parameter) instead 
of a ^/iree- dimensional one. All this simplifies both the 
theoretical treatment and the numerical simulation. For 
clarity, we will mostly discuss our results in the rescaled 
variables (fT2]) - and only occasionally return to the origi- 
nal variables, mostly for comparison with the experiment 
and when we take into account anisotropy. 

In Theorem 3 in jCAOSOTl . p.233], we rigorously show 
that the reduced energy functional is the scaling limit 
of the renormalized full micromagnetic energy in Regime 
III. 



I. Bifurcation 

We now return to the issue of the type of bifurcation 
on the level of the reduced model. Let us note that the 
Hessian of the reduced model in 7712 = can be explicitly 
diagonalized and the first unstable mode is given by in^ = 
sin(7rx2) sin(27r^), where w'' = (327r)-^/^ in agreement 
with (|6]). The reduced critical field is given by 



h* 

'^ext 



3(1) 



4/3 



(13) 



In order to determine the type of bifurcation, we have 
to investigate the energy functional E close to the one- 
dimensional subspace {A1712} generated by the unstable 
mode 771^ . Because of the invariance of both the energy 
E and the unstable mode 77^2 under the transformation 
(7712 —7712 and X2 ^ 1 — X2), all odd terms in the 
amplitude A in the expansion of E{Arh2) vanish. In par- 
ticular the cubic term vanishes so that the bifurcation is 
degenerate. 

This degeneracy of the bifurcation means that at the 
critical field strength /ijxf ^^^^ non- vanishing term 
in the expansion of £^(^4771^) with respect to A is at least 
quartic. Hence it is not sufficient to consider E just along 
the linear space {A7712} but it has to be analyzed along 
a curve {A7712 + ^^7712*} in configuration space. Indeed, 
the curvature direction 77^2* affects the quartic term in 
the expansion and has to be determined such that E is 
minimal. This minimization problem of the coefficient 
of the quartic term is quadratic in 77^2* and thus can be 
solved explicitly. We obtain 

= -ii;(|)'/'sin(27rx2)sin(47rf^), 

which leads to a negative coefficient of the quartic term 
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in the expansion of namely 



(14) 

The negative quartic coefficient implies that the bifurca- 
tion is subcritical or of first order. Subcriticality entails 
that close to m2 = 0, there are only unstable station- 
ary points for he^t slightly below /i*^^, and no stationary 
points close to m2 = for /lext slightly above h^^^^ cf. 
Figure m 



energy 



r 



< K 



\ configuration 
space 



he^t > K 



FIG. 12: Energy landscape close to the bifurcation. The loss 
of stability at the critical field leads to a first-order phase 
transition - on a large scale however the energy is coercive. 



the reduced energy is coercive for all valu es of the exter- 
nal field /lext? see Theorem 4 in jCAOSQTl . p. 236]. This in 
particular implies that there always exists a global min- 
imizer of the reduced energy - which corresponds to a 
local minimizer of the original energy ([T]), see Theorem 5 
in |CA0SQ7l . p. 237] - in particular for fields larger than 
the critical field. But it is not immediately clear how 
these minimizers relate to the unstable mode. 

It is natural to resort to numerical simulations; de- 
tails on the discretization and the algorithms are pro- 
vided in Section [iVl To confirm the conjecture that the 
unstable mode in Regime III is indeed related to the 
concertina pattern, we use a numerical path-following 
algorithm in order to compute the bifurcation branch. 
Figure [14] shows the outcome of the numerical simula- 
tions. As expected due to the coercivity of the energy 
functional, we find a turning point as we follow the bi- 
furcation branch. The turning point is located at a field 
which is just slightly - about one percent - smaller than 
the critical field. After the turning point the configura- 
tions become stable, at least under perturbations of the 
same period. 



At first sight it is surprising that the stray-field energy 
contribution to which gives rise to the only quartic 
term in m2, and clearly is non- negative, may nevertheless 
allow for a negative coefficient in front of the quartic term 
in the expansion ([T4|) . This comes from the fact that the 

two terms in the charge density — + 92m2 interact, 
giving rise to a cubic term in m2 (quartic in A), which in- 
deed allows for cancellations. The way how this operates 
is better understood in physical space: The term mj* in 
+ A^mJ* (the curvature direction in configuration 
space) induces a tilt of the symmetric charge distribution 
of , see Figure [131 This tilt brings opposite charges 
closer together, thereby reducing the stray field energy - 
while increasing the exchange energy to a lesser amount. 






FIG. 13: Unstable mode {A?fi2} and additional curvature 
correction {Arri2 + A^?fi2*} with generated charges. 



FIG. 14: Numerical simulations: The i()*-periodic branch 
close to the bifurcation and the pattern at the indicated fields. 
The gray scales encode the 7712 component but are not compa- 
rable. The whole spectrum is exhausted so that the structure 
of the pattern can best be resolved. By (m|)"'^/^ we denote 
the spatial root mean square of ?fi2, i.e., the amplitude of the 
average magnetization. 



Since the bifurcation is first order, it is not obvious 
whether minimizers of the reduced energy functional can 
be related to the unstable mode. In particular, this find- 
ing sheds doubt on our hypothesis that the concertina 
pattern inherits the period of the unstable mode. It is 
even not obvious whether minimizers of the reduced en- 
ergy functional exist at all. However, one can show that 



As the field increases beyond the turning point, the 
unstable mode develops into a domain pattern of con- 
certina type with its typical scale separation between the 
wall width and the domain size, cf. Figure [151 We thus 
find a continuous transition from the unstable mode to 
the concertina pattern - confirming our hypothesis. 

The numerical simulations lead to the conjecture that 
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FIG. 15: Numerical simulations: The -periodic concertina 
pattern exhibits a clear scale separation as hext increases; 
he^t = 6.84, 23.2, 40.1, 57.3 from left to right. The gray scales 
linearly encode the m2-component and are comparable. 

the magnetization in a perfectly homogeneous, isotropic 
sample exhibits a first-order phase transition from the 
uniformly magnetized state to the concertina state of 
period k;* at the critical field. Clearly this does not 
explain the deviation of the average wavelength in the 
experimental measurements from the theoretical predic- 
tion. Before we address this deviation we now introduce 
a sharp-interface model, so-called domain theory^ that 
is used to investigate the further transformation of the 
concertina for large external fields /lext ^ 1, in particu- 
lar the coarsening, see Subsections IIII A| IIII B\ and IIII Cl 
Our explanation of the coarsening will also provide an 
understanding of the initial deviation of the period, cf. 
Subsection [mDl 

II. DOMAIN THEORY 

In the numerical simulations, we observe for large ex- 
ternal fields a clear scale separation between domains, 
where the magnetization is almost constant, and walls, 
in which the magnetization quickly turns, cf. Figure [151 
This suggests the application of a sharp-interface model, 
namely domain theory. In the following we first discuss 
admissible Ansatz functions and then derive their energy 
within domain theory. This leads to a model which only 
depends on a small number of parameters in configura- 
tion space that is used in Section [iTT] and (Vl] in order 
to get a better understanding of the coarsening of the 
concertina. 
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FIG. 16: Domain theory: The mesoscopic charge-free Ansatz 
function. 

On a mesoscopic scale, the computed magnetization 
is close to a piecewise constant magnetization of ampli- 
tude 7712, ^2 = ±7712 quadrangular domains 
and 7712 = in the triangular domains as indicated in 



Figure [161 We observe that the angles in the pattern are 
related to the amplitude of the magnetization 7772, cf. Fig- 
ure [151 approximately we have that sin a = 2771^. This IS 
related to the fact that the (reduced) stray-field energy 
is strongly penalized for large fields, as we shall explain 
now. In fact, the piecewise constant magnetization is a 
distributional solution of 

-Si(^)+S277l2=0, (15) 

which means that the normal component of the vector 

field (—^,7712) is continuous across the interfaces - this 
is a version of the Rankine-Hugoniot condition in the 
theory of conservation laws. This condition obviously 
holds in case of the vertical walls. In case of the diagonal 
walls the condition 

0=[P.{-^, m2)] = V ■ {-\{ml)\rn%\ (16) 

where v denotes the normal of the diagonal wall as de- 
picted in Figure [TH is equivalent to sin a = 27712- There- 
fore the piecewise constant magnetization satisfying ([16]) 
mesoscopically carries no stray-field energy. Of course, 
on a microscopic scale equation ([T5]) does not hold: The 
continuous transition in the wall generates a right hand 
side, i.e., dipolar charges. Note that walls have to form 
since ([T5]) does not allow for non-trivial smooth solutions 
with boundary data 7712 = 0. 

Within domain theory we therefore consider piecewise 
constant magnetizations of concertina type of period w 
and of amplitude 7772 = ±7713 in the quadrangular and 
7712 = in the triangular domains, s.t. ([16]) holds. Since 
the angles are fixed by ([16]) , admissible configurations are 
characterized by two parameters, namely the amplitude 
of the magnetization 77I2 and the width of the folds w. 
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FIG. 17: Domain theory and numerical simulations: The do- 
main theoretic prediction for the optimal amplitude (red) and 
the computed amplitude based on the reduced model (blue). 
For the reduced model we display the amplitude, i.e., the 
maximal value which is attained in the quadrangular domain. 

The energy which discriminates between these solu- 
tions is given by the total wall energy, which is an appro- 
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priate line-energy density integrated over the interfaces, 
augmented by Zeeman energy. The specific line energy is 
a function of the jump [1712] = of the magnetization 
across the wall - an infinitesimal version of the wall an- 
gle. Due to the shear invar iance of the reduced energy, 
namely 

f 1 = 8X2 + ^1, X2 = X2, 7712 = 1712 - (17) 

by which a diagonal wall can be transformed into a ver- 
tical wall - for the choice of 5 = ±^ - the specific line 
energy can be obtained by restricting the reduced energy 
functional to one-dimensional transitions with boundary 

data ±7722 ^^^^ ^f the vertical walls and ±^ in case 
of the diagonal walls. The optimal transition layers are 
low-angle Neel walls whose line-energy density scales as 

ewalU^j ewall(^2) ^ f (^2) 1^ ^ 

^core 

where Wtaii and Wcore a re the two characteristic length 
scales of the Neel wall, jDKM0Q5L Section 6]. The tails 
of the Neel wall decay only logarithmicly and spread as 
much as possible. In case of the concertina pattern, 
they are only limited by the neighboring walls - thus 
'^taii ^ J . A more careful inspection shows that the core 
width decreases with increasing jump size, more precisely 
'^core ^ (^2)~^7 |SteQ6L Subscctiou 3.5.5]. Hence we 
obtain 

ewau(mO) « ^{ml)Hn-\w{mlf). (18) 

With the rescaling (|9|) and (fTQ|) this turns into 

ewaiiK) « lf{m''^)Un-\d-Hw{m°^f). (19) 

Within the class of admissible magnetizations, the do- 
main theoretic energy becomes a function of only three 
parameters, namely w and /lext- To see this, notice 
that one period of the pattern in Figure [16] contains 

• two vertical walls of height 1 — ^ and of jump size 

7712 

2 7712, leading to an energy contribution of 2(1 — 

;^)ewaii(m^), 

• four diagonal walls of projected height and of 

7712 

jump size 7^2, leading to an energy contribution of 
4 4|e.an(f), 

- 2 

• two quadrangular domains of total area w — 
leading to a Zeeman energy of — /igxt (^2 ) ^ ('^ ~ ^ ) • 

Hence, the total domain energy per period in rescaled 
variables is given by 

^domain(^2 5^ext, U)) 

= 2(1-^) ewaii(m^) + 4 ^ e^aii (#) 
-k^t{ml)\w-^J. (20) 



Within the original scaling the domain theoretic energy 
takes the form of 

-^domain (7712, /^ext,^) 

= 2(£-^)ewan(m^) + 4^ewan(^) 
-he^t{m'2?t(wi-^y (21) 

First of all we apply ([2Q|) to derive the optimal amplitude 
of the i()*-periodic concertina pattern as a function of 
the external field /lext by optimizing the energy in 7712- 
Of course, domain theory is only applicable and thus a 
good approximation for the reduced model for he^t ^ 1 
in which case there is a clear scale separation between 
walls and domains. Figure [171 shows that in this case 
domain theory is in good agreement with our numerical 
simulations. 

Before we go on with the analysis of domain theory 
let us emphasize that the experimentally observed con- 
certina is of course not of uniform period and equal am- 
plitude as our domain theoretic Ansatz above. As shown 
in Figure [TH there are also oblique piecewise constant 
weak solutions of ([T5]) . Nevertheless this class of Ansatz 
functions is very rigid: An elementary calculation shows 
that the location of the interior triplet Aq is uniquely 
determined by the jump condition ([T6|) . if the distance 
between the boundary triplet A2 and Ai, and ml and 
7712 both sides are given. Hence the continuation of 
the pattern is uniquely determined if either the ampli- 
tude in the next quadrangular domain or the location of 
the next triplet, i.e., the width of the next quadrangular 
domain, is prescribed. 



A3 




FIG. 18: Domain theory: Generalized tilted Ansatz function. 



III. COARSENING OF THE CONCERTINA 
PATTERN 

A. Domain theory: The optimal period of the 
concertina pattern for large external fields 

Experiments show an increase in the average con- 
certina period w as the external field he^t is increased 
after the pattern has formed, see Figure [2 The gen- 
eral tendency that the optimal period w is an increasing 
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function of /lext can be understood on the basis of do- 
main theory in the reduced variables he^t and w. By 
optimizing the energy per unit length with respect to the 
period w and the amphtude m2 of the transversal com- 
ponent, we obtain the following scaling of the optimal 
period of the pattern as a function of the external field 

^a(^ext) ~ ^ext In /igxt ^ext > 1- (22) 

In particular we find that the optimal period increases 
with increasing field /lext - the a in Wa stands for absolute 
minimizer. Domain theory also yields the (same) scaling 
behavior for the optimal transversal component of the 
magnetization 

^2a(^ext) ~ ^ext In /igxt ^ext > 1- (23) 

We note that both scalings have also been confirmed by a 
rigorous asymptotic analysis of the reduced energy func- 
tional (fT2]) which does no t rely on a simple concertina 
Ansatz, cf. Theorem 1 in jOSlOL p. 147]. Moreover, nu- 
merical simulations of the reduced energy show that the 
optimal period increases with /lext also for external fields 
close to the critical field, see Figure [191 The optimal 
period shown in this diagram was computed by minimiz- 
ing the energy per unit length both with respect to the 
magnetization and the period, for varying external field. 
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FIG. 19: Numerical simulations: The optimal period of the 
concertina pattern as a function of the external field computed 
on the basis of the reduced model. 



B. Coarsening: A modulation instability 

Although the above analysis predicts that the optimal 
period Wa increases as the field /lext increases, it does not 
explain why and in which way a concertina pattern of pe- 
riod w becomes unstable as /lext increases. We will see 
that both the increasing period for large fields and the 
deviation of the initial period close to critical field from 



the one of the unstable mode are due to an instability 
under long wavelength modulations of the pattern. The 
mechanism behind the instability is the following: Given 
/lext and a period an optimization in the transversal 
component m2 yields that the optimal energy per pe- 
riod Eopt{he^t^w) is a concave function in w if /igxt is 
sufficiently large. The concavity suggests - as depicted 
in Figure [20] - that the concertina pattern of a uniform 
period w becomes unstable towards perturbations which 
increase the period to w-\-e and the corresponding ampli- 
tude of the transversal component to m2{w + e) in some 
folds, and decrease the period to w — e and the amplitude 
to 7fi2{w — e) it in other folds. This modulation eventually 
leads to the collapse of the smaller folds, i.e., the coars- 
ening. However, in view of the non-local character of the 
stray-field energy, it is not clear whether this simplified 
picture, i.e., that the energy of the modulation amounts 
to the modulation of the energy, applies. As we shall see 
in Subsection nil G[ a modulation of the period on a very 
long length scale overcomes this objection. Thus the con- 
cavity of the minimal energy implies an instability under 
long wavelength modulations of the pattern. 




FIG. 20: Concavity of the minimal energy per period implies 
an instability under modulation of the wavelength. 

In order to derive the concavity of the minimal energy 
we apply domain theory for large external fields in Sub- 
section nil CI and an extended bifurcation analysis close 
to the critical field in Subsection IIIIDl We will see that 
both asymptotics match the results of the numerical sim- 
ulation of our reduced model. 

Let us mention that the modulation instability of the 
concertina pattern is closely related to the so-called Eck- 
haus instability which was discovered in the context of 
non-linear instabilities in convective systems leading to 
a change in wavelength of the observed periodic pattern, 
cf. [Eck92j . 

C. Bloch-wave theory: Instability with increasing 
field 

As indicated above, not only the optimal period 
but also the coarsening can be explained on the ba- 
sis of domain theory for large external fields /lext ^ 
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1. This relies on the optimal energy per period 
min^2 ^domain (^2, ^ext, u))- For periods w much smaller 
than the optimal period at some value of the external 
field /lext, i-e., w <C /lext In /lext , we find that 

mini:domain(^2, /^ext, w) -hl^^w'^ ln(/lext^^). (24) 

rh2 

In particular, the optimal energy per period in is 
concave in the period w. Although domain theory there- 
fore suggests an instability under wavelength modulation 
for periods which are much smaller than the optimal pe- 
riod, it is too rigid to allow for such a type of perturba- 
tion, even in the class of generalized Ansatz functions, cf. 
Figure [m 

It is rather on the level of the reduced model that it 
can be seen that the concavity translates into an insta- 
bility (despite the potentially long-range interactions via 
the stray field). Indeed, a so-called Bloch-wave analysis 
of the reduced model shows that the concavity is in a one- 
to-one correspondence with an instablitity under long 
wavelength modulations of the pattern. In the Bloch- 
wave analysis one considers 7Vi()-periodic eigenfunctions 
of the Hessian of the form 

57712 = e-'^'^' 8ml' 

with wavenumber ki = and N some large integer 

and where Srril^ is i()-periodic with respect to fi, i.e., 
one condisders sinusiodal modulations of some suitable 
t()-periodic function. An asymptotic expansion of 

Hessi:(m2)(e-^^^^Mm^^) = A^^e"^^^^^ 5ml' (25) 

for small wavenumbers /ci <C 1, i.e., TV ^ 1, shows that 
the first eigenvalue can be related to the second derivative 
of the optimal energy per period £^opt = min^2 More 
precisely, one can show that the eigenvalue possesses the 
expansion 

. d^ ^ ^ 
^ Co kl-^Eo^t{h^^t,yo) for A;i <Cl, 

where Co denotes a constant that depends on m2, see 
[StelOl Theorem 5.1]. This shows that the concavity of 
Eopt{hext^w) with respect to the period w implies that 
the concertina pattern of a given period w is unstable. 
Domain theory predicts that the marginally stable pe- 
riod Ws, i.e., Ws such that ^£^opt(^ext, ^s) = 0, scales 
as Ws ~ /^ext In ^ext 7 cf . ([24|) - we note that the s in Ws 
stands for marginally stable. Figure [2T] displays the op- 
timal and the marginally stable period computed on the 
basis of the reduced energy functional. 

Figure [22] shows that the computation of the optimal 
and the marginally stable period on the basis of domain 
theory matches the numerical simulations on the basis of 
the reduced model. 




FIG. 21: Numerical simulations: Comparison of the optimal 
and marginally stable period of the concertina pattern as a 
function of the external field - both computed on the basis 
of the reduced model. In the region below the red curve the 
minimal energy per period is concave and thus a concertina 
of that period is unstable and coarsens. 
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FIG. 22: Numerical simulations and domain theory: The op- 
timal and marginally stable period computed on the basis 
of the reduced model (dashed) match the predictions on the 
basis of domain theory in the regime /lext ^ 1- 



D. Bifurcation analysis: Instability for small fields 

The numerical computations, cf. Figure [STJ show 
that the optimal energy per period is concave not only 
for large external fields as predicted by domain the- 
ory. In fact, we extract from our numerical data that 
^2^opt(^ext7 '^*) is negative also for small external fields 
up to the turning point. This is consistent with the 
numerical computation of the eigenvalue based on 
the asymptotic expansion of equation ([25|) . Hence, the 
Bloch-wave analysis implies that the -periodic con- 
certina pattern is unstable under long wavelength mod- 
ulations close to the critical field. 
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This qualitatively explains the trend in the deviation 
of the initial concertina period w ^^^ f rom the period of 
the unstable mode, see Subsection II Gl Close to the crit- 
ical field, the concavity can be confirmed with the help 
of an asymptotic bifurcation analysis. To see this, we 
extend our Ansatz from Section [Hi and take into account 
small deviations of the wavenumber k = -\- Sk. As we 
have seen in ([T4|) in Section [Hi the quartic coefficient in 
the energy expansion, namely g^^, is small compared to 
the second order coefficient and the scale of the reduced 
external ffeld. Due to that degeneracy it is necessary 
to additionally take into account a contribution of cu- 
bic order in the perturbation of m2 = 0, i.e., we use the 
extended Ansatz 

7712 ~ + A^m2* + A^mJ**. 

Optimizing the coefficients in the expansion of the energy 
with respect to A subsequently in mj* and m2** leads to 
an expansion of the energy density of the form 

^i(^ext(^) - he.t)A^ - c^{k)A^ + C6(^)A^ 

where C4(/c*) = in accordance with (fH|) . Hence 

under the assumption that C4(/c*) ^ 0.00105 is small, the 
energy density to leading order can be approximated by 

- + + c,{t)A'. (26) 

The numerical values of the coefficients are given by 

i^4(^)|^.^*^ -0.0217, 
C6{t) ^ 0.000207. 

Notice that C6(/c*) is positive, confirming the numerically 
observed turning point of the i()*-periodic branch. Obvi- 
ously, the asymptotic expansion displays an asymmetric 
behavior in Sk; the energy decreases for Sk < 0. Based 
on the expansion ([26|) . one can characterize the optimal 
wavenumber and the optimal period. We note that the 
concavity of the minimal energy per period as a function 
of the period is equivalent to the concavity of the energy 
density as a function of the wavenumber k: 

^^(-) = (&S(^^(t))- 

Figure [23] shows the optimal period and the marginally 
stable period calculated on the basis of ([26]) . We read off 
that the i^)* -periodic concertina pattern is indeed unsta- 
ble at the critical field. 




FIG. 23: Bifurcation analysis: The optimal and marginally 
stable period as a function of the external field obtained on 
the basis on the extended bifurcation analysis. 

A comparison between Figure [2T] and Figure [23] shows 
that the predictions on the basis of the asymptotic ex- 
pansion differ from the optimal and the marginally sta- 
ble period computed on the basis of the reduced model, 
compare for example the scale of the external field. This 
deviation is related to our assumption that the quartic co- 
efficient is small so that the energy can be approximated 
by ([26|) . On the other hand. Figure [24] shows that the 
asymptotics match the reduced model if we add a quar- 
tic contribution +■§ / to the reduced energy where 
the value of the parameter Q is such that the contribu- 
tion cancels C4(^*) in ([26|) ( which happens Q ^ 0.03). 
We will see later that such an additional quartic contri- 
bution has a physical meaning if we take into account a 
uniaxial anisotropy, see Section [Vl] It turns out that Q 
corresponds to an appropriately rescaled quality factor 
Q. 



E. Numerical bifurcation analysis: Type of 
secondary instability and downhill path in energy 
landscape 

With the help of a bifurcation-detection algorithm we 
are able to compute at which field the i()*-periodic con- 
certina becomes unstable under A/" iD* -periodic perturba- 
tions while we follow the primary branch. Figure [25] 
shows the secondary critical fields; as expected (cf. Sub- 
section [TTlD] and Figure [2T]) the secondary instability ap- 
proaches the turning point as the integer N increases. We 
note that it is reached for finite N. 

In the following we want to study in which way the con- 
certina pattern becomes unstable. We first present the 
outcome of the computation of the secondary bifurcation 
branches. We point out that due to the symmetries of 
the pattern, the bifurcations are not simple in the sense 
that more than one branch bifurcates. 
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FIG. 24: Numerical simulations and bifurcation analysis: 
The prediction on the basis of the reduced model (dashed) 
matches the prediction on the basis of the extended bifurca- 
tion analysis for a near- degenerate value of Q = 0.0295 close 
to Q* ^ 0.03, cf. Section Ell 
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FIG. 26: Numerical simulations: Bifurcation diagram for 2w* 
perturbations: The bifurcation branches that connect the w*- 
periodic (blue) and the 2it)* -periodic branch (orange). The 
magnetization patterns at the indicated fields are shown in 
Figure [271 
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leads to the collapse of exactly one fold - reducing the 
total number of folds from TV to — 1, see Figure [25l 
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FIG. 25: Numerical simulations: The appearance of the sec- 
ondary instability under A/'if)* -periodic perturbations as a 
function of N. The critical field for A/" = 8 is given by 5.602. 

The symmetries of the pattern can be identified as lin- 
ear representations of the dihedral group 1^2 at, where TV 
indicates the number of folds. The secondary bifurca- 
tion branches are computed with the help of a numerical 
branch switching algorithm which is adapted to the prob- 
lem of multiple bifurcations. Generically, there are two 
distinct types of branches: Branches along which rota- 
tional symmetry is broken and reflectional symmetry is 
conserved and vice versa, see Figure [27l In case of the 
first type of branches, a fold collapses as two neighboring 
faces disappear; in case of the second type of branches, 
the number of folds decreases as one face disappears and 
the two adjacent faces merge. During the coarsening pro- 
cess, the width of the remaining folds is adjusted. Let 
us point out that the first instability of the t()-periodic 
concertina under A^i()-periodic perturbations in the end 



FIG. 27: Numerical simulations: Reflectional symmetric with 
respect to center wall (top) and rotational symmetric with 
respect to the midpoint of white face (bottom) magnetization 
pattern on the unstable bifurcation branch connecting the 
If)* -periodic and the -periodic branch. The central fold 
collapses (top); white face disappears and two adjacent black 
faces merge (bottom). 



F. Wavelength modulation in the experiments 

In the experiments, the xi-wavelength of the modu- 
lation is restricted by the finite extension of the sam- 
ple. Moreover, inhomogeneities and defects of the mate- 
rial, in particular those at the edges of the cross section, 
strongly affect the formation. This is reflected by the 
fact that walls occur at the same pinning sites when the 
experiment is rerun. The existence of pinning sites hence 
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FIG. 28: Experiment and numerical simulations: The coarsening of the concertina pattern in a Permalloy sample (top row) of 
30nm thickness and 70|am width compared to the numerical simulations (bottom row). A ripple- like structure grows into the 
concertina pattern. Within the numerical simulations we iteratively increment the external field and minimize the energy. The 
computational domain is of period 6w* . The numerical images are scaled according to (|9]). The numerical image hence display 
1.8 times the unit cell; the numerical images therefore appear to be more uniform than the experimental concertina. 



leads to an effective modulation wavelength that is just 
a small multiple N of the wavelength of the pattern. In 
particular we expect that pinning sites have a stabilizing 
effect and therefore prevent coarsening. Therefore, the 
seemingly artificial numerical simulation for small and 
moderate cf. Figure [25l may be more relevant for the 
experiment than the Bloch-wave analysis, i.e., N oo^ 
cf. Section [mCl 



G. Domain Theory: Instability for decreasing field 

The experiments also show that the concertina period 
w decreases with decreasing external field /lext- This 
has a simple explanation on the level of domain theory, 
too. Suppose that the concertina period had increased 
at several coarsening events during the increase of the 
field. As the decreasing external field /igxt drops be- 
low its optimal scaling given the period that is, for 
w ^ /lext In ^ext7 the optimal concertina pattern does not 
suffer a long wavelength instability, but instead degener- 
ates in the sense that the closure domains invade the 
whole cross section. Simulations of the reduced model 
confirm this scenario predicted by domain theory, see Fig- 
ure [29l which shows a pattern of period bw"^ close to the 
turning point: The numerical backward cycle, in which 
we start at the multiply coarsened state and then after 
minimization repeatedly decrease the external field by a 
fixed increment, shows that the coarsened pattern stays 
stable up to the turning point that coincides with the 



moment at which the pattern degenerates as mentioned 
above. Depending on the initial level of coarsening, the 
period is then either refined or we reach the uniformly 
magnetized state after the minimization. 
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FIG. 29: Numerical simulations: The coarsened concertina 
pattern degenerates as the external field is reduced. The nu- 
merical simulations confirm the prediction based on domain 
theory: The pattern degenerates at the turning point of the 
branch. 



H. Conclusion: Hysteresis and scattering of data 



Summing up, domain theory in conjunction with a 
Bloch-wave argument indicates that the concertina pat- 
tern of period w is present or stable at a given field /lext if 
and only if ^ - £H-^he^tln{d-^/^t/H-^^^he^t), which 
is confirmed by the numerical simulations. In particu- 
lar we expect that the height of the triangular domains 
(~ ^) is close to constant as the external field increases, 

cf. ([2^ and ([23]) . If the period deviates by a (large) 
factor from that expression, it becomes unstable. On 
the other hand, this analysis also suggest that there is 
a range of w ^ £^t~-^/iext In (i~^/^-^^/^t~^/^/iext for which 
the concertina pattern is stable, see Figure [3Ql This may 
explain some of the scatter in the experimental data and 
the pattern's hysteresis. 
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FIG. 31: Numerical simulations: The hysteresis loop. As we 
increase the external field we follow the red path: The con- 
certina pattern coarsens if the period is smaller than the stable 
period. As we decrease the field we follow the yellow path: 
Starting from a coarsened concertina the pattern degenerates 
as we reach the turning point of the branch. The pattern 
refines towards the optimal period until it finally disappears. 




FIG. 30: Numerical simulations: The marginally stable (red), 
optimal (blue) and maximal period (green) of the concertina 
pattern as a function of the external field /lext- The dashed 
and solid curves depict the result on the basis of the reduced 
model and on the basis of domain theory, respectively. 



Figure [30] displays the marginally stable period (below 
the red curve the minimal energy per period is concave 
and thus the concertina of smaller period unstable as the 
field increases) and the optimal period depending on the 
external field. The upper green curve indicates the turn- 
ing points of the i^)m-periodic branches, i.e., the smallest 
external field for which a concertina of a certain maxi- 
mal period Wm exists - clearly the m in Wm stands for 
maximal. Observe that the maximal period Wm on the 
basis of domain theory and on the basis of the reduced 
model coincide for large external field, too. The region 
bounded by Wg and Wm corresponds to the range of stable 
periods. 



IV. DISCRETIZATION AND NUMERICAL 
SIMULATIONS 



The numerical simulations are based on a finite differ- 
ence discretization of the reduced rescaled energy func- 
tional ([T2]) . The transversal component m2 is approx- 
imated on a uniform Cartesian grid. The discretiza- 
tion of the exchange, anisotropy and Zeeman energy is 
straight forward. In case of the non-linear charge density 

a = —di^-\-d2rh2 our choice of a finite difference sten- 
cil is motivated by the inheritance of the shear invariance 
([TT]) . The stray-field energy can efficiently be computed 
using Fast Fourier Transform with respect to xi. For 



an introduction of the discretization scheme see iSteOGl . 
Subsection 3.2]. Note that the computation of the energy 
and related quantities, such as gradient or Hessian, can 
be parallelized since the non-locality is only with respect 
to one dimension. For the parallelization we decompose 
the computational domain into horizontal slices with re- 
spect to X2. 

We apply numerical simulations to compute (local) 
minimizers and stationary points. The naive approach 
using steepest descent algorithms for the computation 
of minimizers is slow and even fails close to bifurcation 
points. The iterative path- following techniques that we 
apply in order to compute an approximation to a branch 
of stati onary points are adapted to such situations, cf. 
|GeoOl( . The local tangent in a stationary point 
(m^, /^ext) ^f branch is used to obtain a predictor for 
the next point on the branch (^2^^, /i-ext^)? ^f. Figure 
[33l Within the corrector step the predictor is orthogo- 
nally (to the tangent) projected onto the branch. This 
step amounts to the solution of a non-linear equation. 



18 




FIG. 32: Experiment: The hysteresis cycles of a Permalloy sample of 30nm thickness and 50|J,m width. The upper row shows 
the pattern as the external field increases (from left to right), the lower row shows the pattern as the external field decreases 
(from right to left). 
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FIG. 33: Tangent predictor-corrector continuation method. 



more precisely an augmented Euler-Lagrange equation: 



= 0. 



A bifurcation point can be detecte d with the help of 
an appropriate indicator function, cf. |GeoQl| . However, 
both the bifurcation detection and the branch-switching 
technique which are described in that reference are appli- 
cable for simple bifurcations points only. As described in 
detail in both methods can be modified in order 

to cope with multiple bifurcation points. This extension 
relies on the fact that multiple bifurcations which occur 
due to symmetries of the primary solution generically can 
be reduced to simple bifurcation points, cf. [GSQ2i] . 



V. POLYCRYSTALLINE ANISOTROPY 

The experiments usually do not show a clear-cut crit- 
ical field with a first-order transition (i.e., subcritical 
bifurcation). This can be due to lack of experimental 
resolution (the amplitude of the transversal component 
(m|)^/2 = (p/^£-^/H-^^^{miy^'^ at the turning point 
ranges between 0.015 to 0.063 for typical sample dimen- 
sions, namely widths i between 10 to 50|LLm and thick- 
nesses t between 30 and lOOnm) or due to the presence of 
the so-called ripple that smoothes out the transition, as 
we shall explain in this section. The ripple is the in-plane 
small-scale oscillation of the magnetization - perpendic- 
ular to its average direction - in extended films. In this 
sec tion, we show how the linear ripple theory developed 
m jHof68LlHar68l can be incorporated into our theory for 
the concertina - and explains the smoothing-out of the 
first-order transition encountered in Section [iTl 

The ripple is triggered by an effective field of random 
direction on a small scale. Several origins for this effective 



field ha ve been proposed in the literature, see for instance 
|Har68l Section C]; in polycrystalline thin films, the ran- 
dom orientation of the grains (via crystalline anisotropy) 
and local stresses (via magnetostriction) are seen as the 
main causes. In our discussion, we focus on the former. 

Hoffmann |Hof68 | and Harte |Har68j . based on the 
torque equilibrium, linearized around a spatially con- 
stant magnetization (solely determined by the external 
field and anisotropy). Hereby they identified the lin- 
ear response to (for instance) such a small-scale, small- 
amplitude random effective field. The main finding is 
that the stray field - which penalizes transversal more 
than longitudinal perturbations of the magnetization be- 
cause the former lead to a stronger charge oscillation 
- resu lts in a strong anisotropy of the response. Hoff- 
mann |Hof68 | characterized this respons e in terms of the 
Green's function, whereas Harte |Har68( characterized it 
in terms of the multiplier in Fourier space, i.e., /c-space. 
Since Hoffmann chose to expand the Green's function (in 
a self-consistent way on the leve l of th e length scale) in 
terms of a Bessel function, see |Hof68l (5)], it devia tes 
order one from the exact expression in jKre67LlBro70j . 

Clearly, the anisotropic rescaling (j9|) leading to our re- 
duced model and the anisotropic response have the same 
origin. We will see that both the ripple and the transition 
between ripple and concertina can be explained within 
the framework of an extension of our reduced model. We 
note that our analysis of the ripple is mainly a reformu- 
lation of the classical results. However, the new insight 
is that the finite width £ of the sample leads to a (con- 
tinuous) transition from the ripple to the concertina. 

We now explain how to extend our reduced model. We 
start from the 3-d model ([T]) with a uniaxial anisotropy 
of strength Q and position-dependent easy axis e(x), i.e., 
with the term —Q J{m ■ e)"^ dx. In the approximation of 
our reduced model, i.e., ma = 0, m = m(xi,X2), and 

2 

the linearization mi ^ 1 — ^ due to m2 <C 1, this term 
is, up to additive constants, to leading order approxi- 
mated by —2Qt J m2eie2 dxi dx2, where 6162(^1,^2) de- 
notes the vertical average of the product of the first two 
components of the easy axis e = (61,62,63). A random 
anisotropy therefore acts to leading order as a random 
transversal external field 



2t 



JVL' 



pple^2 dXidX2, 



(27) 
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where /irippie = QeTe2- As mentioned, the position de- 
pendence of e arises from the random orientation of the 
grains of size i grain- Provided t <C i grain <C w"^ (where 
we take as a typical length scale of the magnetization 
pattern), the stationary statistics of 6162 are character- 
ized by 

(eIe^(0,0)eTe^(xi,X2)) 

= Cam^(^i)^(^2)(eTe5(0,0)2), (28) 

where (•) denotes the ensemble average and S the Dirac 
function. 

For subcritical fields /igxt < ^2xt5 neglect the non- 
linear term in the stray-field energy in ([7j). The resulting 
energy functional is quadratic and linear in m2, hence it 
is conveniently expressed in terms of ^m2(/ci, /C2), which 
denotes the Fourier transform of 777-2 in ^1 and the Fourier 
sine series in X2' 
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The explicit minimization yields 

(29) 

We interpret this m2 as the ripple. Since ([28]) on the 
level of TeTe^ reads {\Te]e^{ki,k2)\'^) = £l^^.^, ([29]) is 
best expressed in terms of the energy spectrum: 

(|^.^.(fc,,fc.)P) = Q^ ^^^^^^,J;;^_^^^^^, (30) 

This formula clearly displays the afore mentioned 
anisotropic response of m2 to the isotropic field /^ripple- 

From formula ([30]) one can infer the predominant 
wavenumber of the ripple, that is. 

For moderate stabilizing fields t'^d~'^ —he^t ^ 

^-2/3£4/3^-2/3^ we obtain from §^ that the average 
wavenumber scales as {\ki\) ~ {—hQ^ty^'^d~^ <C td~'^. 
This is the scaling of the pre domina nt wavenumber of 
the ripple in an extended film jHof68L p.34, (7)]. Notice 
that the lower bound characterizing Regime III is equiv- 
alent to t'^d~'^ > For large stabilizing 
fields — /lext ^ t^d~'^ one can show that the average am- 
plitude of the ripple, given by / Xlfe2 ^l-^^^P) d/ci, tends 
to zero. Moreover, from ([3T]) . because of the discreteness 
of /c2, we can infer 

hm (|fei|) = ^, 



which is the wavenumber of the unstable mode ([6]). We 
thus learn that, as the strength /lext of the external field 
increases from negative values towards the critical value, 
the average wavelength of the ripple continuously in- 
creases 

• from the values characteristic for a film which is 
infinite in both xi and X2-directions 

• to the wavelength of the unstable mode that is at 
the origin of the concertina pattern (which depends 
on the sample width). 

Due to this transition it is thus not surprising that rip- 
ple and small-amplitude concertina are difficult to distin- 
guish. 

We now address the numerical simulation of our aug- 
mented model ([32]) . Let us therefore first rewrite the 
additional term ([27]) in the rescaled variables ([9]). The 
rescaled reduced model ([T2]) is augmented by 



ipple"l2 AxiAX2, 



(32) 



where /irippie is a stationary Gaussian field of vanishing 
mean and of variance 

(^ripple(0,0)/lripple(^l,^2)) = (cr* )^(5(f i )(5(f 2 ) , (33) 
With a* = d-^^/H^/H-^/^Q^gra^n{e^{^.^)Y^^• 1^ 

case of a uniform distribution of the anisotropy axis in 
the plane, we have for example that (6162(0, 0)^) = |. 

On the level of the discretization, the field /irippie 
is modeled as a Gaussian random variable of mean 
zero, which is identically and independently distributed 
from grid point to grid point and has variance 
((j*)^Axi~"^Ax2~"^, where Ax^ denotes the grid size in 
direction Xi. For the numerical simulations we thus have 
to determine the value of a* for a typical sample. Let us 
consider a film of 30nm thickness and 70|LLm width with 
typical grain size ^grain = 15nm. For a local strength of 
anisotropy Q = 5 x 10~^ we obtain that (cr*)^ = 125.87. 
For the value of (cr*)^ = 110.83, our numerical simulation 
indeed shows a continuous transition from the ripple to 
the concertina pattern instead of a first-order phase tran- 
sition due to a subcritical bifurcation, see Figure [28l 

We also believe that our reduced model is the appro- 
priate framework to analyze the non-linear corrections 
to the linear ripple theory. Indeed, we have seen in Sec- 
tion [iT] that it captures the transition from the unstable 
mode to low-angle symmetric Neel walls. We thus be- 
lieve it also captures the transition from the ripple to the 
block ed stat e that is related to hysteresis in extended thin 
films jFel6l| . 

Closing this section, we contrast the ripple, that can 
be seen as a consequence of quenched disorder^ to the ef- 
fects of thermal fluctuations. Thermal fiuctuations can be 
modeled by a random external field term in the Landau- 
Lifschitz-Gilbert equation that is white noise in space and 
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time. The reason for modeling thermal fluctuations by 
a space-time white noise torque in the Landau-Lifschitz- 
Gilbert equation is that the stationary measure of this 
Langevin equation is given by the Gibbs distribution 

^exp(-£:(m))dm, (34) 

where E{m) is the 3-d micromagnetic energy functional 
and dm is best thought of as the high, but flnite- 
dimension al mea sure after spatial discretization of m. 
Following jBGQ5L Subsection 2.4], we consider a situa- 
tion where the constant magnetization, say m = (1, 0, 0), 
is a strict global minimizer of E (because of a suffi- 
ciently strong external fleld in direction (1,0,0)). This 
justifles to replace E{m) and dm in (|34|) by its Hessian 
Hess ((5m, (5m) in (1,0,0) and d(5m2d(5m3, respectively. 
This has the advantage that we obtain a Gaussian mea- 
sure that can be expHcitly analyzed. The outcome is 
the following: In case of a bulk material, the expected 
value of + diverges as the mesh size Ax goes to 
zero; the expected value of the wavenumber \k\ behaves 
as Ax~^. The same holds for thin fllms, although the 
divergence is just logarithmic. 

This simple analysis highlights the need of a renormal- 
ization in case of thermal fluctuations. As we have seen, 
quenched disorder coming from polycrystallinity can be 
modeled by a random fleld term in the micromagnetic 
energy that is white noise only in space. As opposed to 
thermal fluctuations, there is no divergence in the ampli- 
tude of the excitations in case of such a fleld term that is 
white noise in space only - the critical dimension for this 
random effect is four. Moreover, in thin films, the dom- 
inant wavelength of the in-plane fiuctuations excited by 
such a field is determined by both exchange and stray- 
field energy and turns out to be much larger than the 
atomistic length scale d and the typical grain size ^grain- 

VI. UNIAXIAL ANISOTROPY 

We now address the effect of uniaxial anisotropy - 
constant throughout the sample - on the formation of 
the concertina pattern. We focus on the two cases in 
which the easy axis coincides with the X2-axis (transver- 
sal anisotropy e = (0,1,0) in ([1])) or in which the easy 
axis coincides with the xi-axis (longitudinal anisotropy 
e = (1,0,0) in ([Tj)). Clearly, such type of anisotropy 
has no effect on the stationary point of the energy, i.e., 
the uniform magnetization. On the level of the reduced 
model both cases can be represented (up to an additive 
constant) by the additional quadratic term 

— Q t j m\ dxi dx2 (35) 

with a signed quality factor Q. Transversal anisotropy 
corresponds to Q > 0, longitudinal anisotropy corre- 
sponds to Q < 0. 



As will become clear below, when considering the ef- 
fects of anisotropy, it is appropriate to expand the Zee- 
man term to quartic order, i.e., 

-he^t t j {ml + ^ ) dxi dx2. 

The following Gedankenexperiment is helpful in under- 
standing the sequel: In extended thin films, i.e., ^ = oo, 
there is no incentive for a spatially varying magnetization 
so that we may consider a constant magnetization m2 in 
which case the relevant energy per volume is given by 

4 

—Qm2 — he^t{m2 + ^). In this case the critical field 
is given by h^^^ = —Q. For longitudinal anisotropy, 
the bifurcation is subcritical, whereas for transversal 
anisotropy, the bifurcation is supercritical and yields 

m2 = ±(2(l + Q-i/iext))'/'. (36) 

Hence for finite ^, there are two competing mechanisms 
which lead to a bifurcation and a selection of an ampli- 
tude for m2 : uniaxial anisotropy and shape anisotropy in 
form of the stray-field energy. 

As we will see in the sequel, there are essentially three 
different effects of anisotropy: a linear one, a weakly 
non-linear one, and a strongly non-linear one, which 
we list and characterize below. Notice that the order 
at which these effects arise with increasing anisotropy 
does not agree with their ordering with increasing non- 
linearity, cf. Figure [34l The linear effect becomes pro- 
nounced for IQI > the strongly non- 
linear one for \Q\ ^ £~^t, and the weakly non-linear 
one only for \Q\ > d~'^^^i~'^^^t^^^ . Note that we have 
that d^/H-^/H^/^ < i-H < d-2/3£-2/3^4/3 provided 
d^l~^ <C t, which is the lower bound on the film thick- 
ness which characterizes Regime HI. 

Strongly non-linear effect 

— \ 1 \Q\ 

^2/3^-4/3^2/3 ^-2/3^-2/3^4/3 

Linear effect Weakly non- linear effect 

FIG. 34: The order of the different effects of anisotropy 

We mainly focus on the case of transversal anisotropy 
Q > 0. In case of longitudinal anisotropy Q < we 
give an explanation for the experimental fact that the 
concertina cannot be observed at all. 

a. Linear effect for weak anisotropy \Q\ ^ 
^2/3^-4/3^2/3 obvious effect of anisotropy is a shift 
of the critical field /i*^^ by the amount — Q; we call it the 
"linear effect" of anisotropy since it arises on the level 
of the linearization at m2 =0. In view of the scaling of 
the critical field /i*^^ at Q = 0, i.e., ([10]), we infer that 
the value of the critical field is dominated by the uniaxial 
anisotropy, i.e., 

K^,^-Q for \Q\^d^^H-^^H^^\ (37) 
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We note that a transversal anisotropy decreases the dis- 
tance between the two critical fields i/i-e^t correspond- 
ing to the stationary states ±m*; in particular, for Q ~ 
^2/3£-4/3^2/3^ the critical field chan ges sign and thus the 
order between the two critical fields switches. (Likewise 
for longitudinal anisotropy the distance decreases.) Al- 
though a clear-cut critical field cannot be observed in the 
experiments due to the polycrystalline structure which 
triggers the ripple, and since the value of the effective ex- 
ternal field at the investigated sample section is not avail- 
able, the linear effect could be qualitatively confirmed: 
For Permalloy samples of high (transversal) anisotropy 
the oscillatory instability occurs before the external field 
is reversed. In accordance with (|37|) . we observe for rela- 
tively wide films that the relative strength of anisotropy 
increases and the critical field decreases (theoretically ap- 
proaching —Q). On the other hand for low-anisotropic 
Permalloy the first oscillation is observed close to zero 
external field. 

h. Weakly non-linear effect for strong anisotropy 
IQI > t{w*)-^ - d-^/^^-^/^t"^/^. For sufficiently strong 
anisotropy Q, the quartic term coming from the stray- 
field energy no longer dominates the quartic term coming 
from the Zeeman energy near the bifurcation. We call 
this effect the "weakly non-linear effect" of anisotropy, 
since it can be analyzed on the level of an expansion 
of the reduced energy near m2 = and /lext = ^Lt^ 
(p!4|) , where we take into account the quartic Zeeman term 
-^tA^ J (m^)^ dxi dx2. The shift of the critical field 
suggests the following rescaling for the reduced external 
field 

/.ext = ^-'/W-2/^(/.ext+g). 

In addition we set 

SO that we obtain with the same rescaling of energy, 
length and magnetization as in (|9|) and (fTTj) the reduced 
energy functional augmented by 

-\-Q / m\dxidx2- 

Therefore the energy close to the bifurcation takes the 
form of 

£:(Am* + A^m**) 

^ -(f )^/'(/^ext - k^,) A' + (^g - ^,)A\ 

For IQI > ^-2/3^-2/3^4/3 > ^2/3^-4/3^2/3 ^^ie Critical 
field asymptotically behaves as h^^^ ^ — Q, cf. (|37|) . 
so that the reduced quality factor behaves as Q ^ 
1 ^2/3^2/3^-4/3 Q ^iQgg ^Yie critical field. From the lat- 
ter we read off that in the regime Q > (i-2/3^-2/3^4/3 
the quartic coefficient becomes positive and therefore the 
bifurcation becomes supercritical, cf. Figure [35l Essen- 
tially it is a perturbation of the constant-magnetization 



bifurcation in infinitely extended films mentioned above, 
cf. (|36|) . In particular, the selected amplitude in this case 
scales as m2 ~ A ~ (1 + /^-extQ"^)^^^- On the level of 
the extended bifurcation analysis one finds that the pe- 
riod of the unstable mode k;* lies in the stable region in 
the neighborhood of the critical field. In agreement with 
this, for increasing external fields the numerical simula- 
tions show that no modulation instability occurs and that 
there is no coarsening. We note that domain theory is 
consistent with the numerical simulations, too. 




^2/3£-4/3^2/3 

FIG. 35: Numerical simulations: Transition from sub- to su- 
percritical bifurcation as strength of transversal anisotropy 
increases. For Q — 0.03 ~ Q* the bifurcation degenerates. 

On the other hand, for large longitudinal anisotropy, 
i.e., -Q > (i"2/3^-2/3^4/3^ expect that there is no 
turning point on the bifurcating branch so that it re- 
mains unstable - to the effect that no concertina pattern 
forms in the first place. The numerical simulations in 
Figure [36] show a second turning point which coincides 
with the break-up of the concertina pattern. For even 
larger longitudinal anisotropy the first turning point is 
destroyed, cf. Figure [36l 

This observation can also be confirmed on the level 
of domain theory where we take into account anisotropy 
and the quartic term in the Zeeman energy, cf. (|2T]) : 

edomain(m^, ^) = 2 - ^) e(m^) + 4 e (^) 

-h,^,\{ml)H(w^-^y (38) 

The quartic wall energy cannot compensate the destabi- 
lizing quartic Zeeman contribution provided he^ttw 

(up to a logarithm). Therefore due to /i*^^ ^ —Q and 
w ~ dP^^P^^t~^^^ close to the bifurcation there are no 
(local) minimizers of the energy. 

Typical values for our Permalloy samples of strong 
uniaxial anisotropy range from Q = . ,_2/3L^2/3^4/3 ^ 
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^2/3^-4/3^2/3 



FIG. 36: Numerical simulations: Loss of the turning point as 
strength of longitudinal anisotropy increases 



2.1 X 10 ^ to 0.023 depending on the sample's width and 
thickness (Q = 5 x 10~^, t = lOnm to 150nm, i = 10|LLm 
to 50|LLm). Typical values for CoFeB range from Q = 
7.8 X xlO-"^ to 0.011 {Q ^ 1.5 x 10"^, t = 30nm-100nm, 
i = 10|LLm-50|LLm). The uniaxial anisotropy is thus too 
small to cause the weakly non-linear effect. However, al- 
though local minimizers of the energy might exist in case 
of longitudinal anisotropy, still the energy is not coercive 
as soon as the external field is reversed. 

c. Strongly non-linear effects for moderate 
anisotropy \Q\ ^ i~^t. In that case one can dis- 
tinguish two different scenarios in the formation of the 
concertina: 

• Scenario I: If the amplitude (and shape) of the con- 
certina pattern would not be affected by anisotropy 
(besides the critical field at which it bifurcates), like 
in an infinitely extended film, its optimal amplitude 
would scale as 



m2a 



1 (|33 

£t-\K^t-K^,) ^ it-\h,^t^Q) 

= ^t-^Q(l + Q-^/iext), (39) 



up to a logarithm for /lext - K^t > d^l '^i-^'H^I'^ , 
as we have seen in ([23]) in Section IIII A[ 

• Scenario II: If the amplitude of the concertina pat- 
tern would be dominated by transversal anisotropy, 
it would behave as 



^2 a 



(1361) 



(1 + g-^/iext)'/' 



(40) 



for 0< (1 + Q-i/iext) <C1. 



Hence we expect that for Q ^ i~^t^ the concertina 
pattern is limited by stray field effects as long as < 
l + (5~^/i^ext <^ {Q~^t~^t^^ and by anisotropy effects once 
(Ci~^t~^tY <C 1 + Q~^he^t < 1- Loosely speaking, the 
effect of anisotropy kicks in for a large amplitude and 



is most prominent close to field strength where the con- 
certina pattern vanishes. We call this the "strongly non- 
linear effect" of anisotropy. (Also this provides a reason 
to expand the Zeeman term to higher order.) 

We note that we have to take into account the lower 
order wall energy in Scenario II in order to determine 
the optimal period. In that case, a minimization of the 
energy per length yields the following scaling behavior of 
the optimal period (up to a logarithm) 



As we know from Section IIIII the experimentally more 
relevant quantity is the marginally stable period, i.e., 
the largest period (as a function of the external field) 
for which the minimal energy is convex. At the cross- 
over we expect that the marginally stable period is 
of the order ~ tQ~^^ cf. Figure [39l In fact, due to 
{Q~^£~^t)'^ ~ 1 + Q~^hext at the cross-over, we have 
that w - ft-^Q{l + Q-^/iext) ^ tQ-^, see §9^} to- 
gether with the fact that Wg ~ a- ^ period 
of that order the minimal energy in Scenario II turns 
out to be convex. Hence we expect that the coarsen- 
ing stops once {Q-^i'^t)'^ <C 1 + Q~^K^t < 1- Still 
the transversal component of the magnetization grows as 
m2 ~ (1 + Q~^hQ^tY^'^ so that size and height of the 
closure domains decrease. 





FIG. 37: Experiment: Permalloy samples of width 60|am to 
150|J,m of high anisotropy and at the end of the coarsening 
process. The 6 samples on the right are of thickness 30nm, 
the 6 samples on the right are of thickness 50nm. The period 
of the pattern appears to be independent of the width of the 
samples, in agreement with our theoretical prediction as an 
effect of anisotropy. 



Figure [39] displays the transition of the scaling behavior 
in the optimal period, the marginally stable period, and 
the amplitude of the transversal magnetization compo- 
nent. At the cross-over we have that m2 ~ t{£Q)~^ <C 1 
and w ^ tQ~^ <^ i. This is consistent with the assump- 
tions of the reduced model, i.e., the low-angle approxi- 
mation and the scale separation of the dominant length 
scales with respect to xi and X2. For the same reason 
and due to the observation that /lext + Q > (f^^i~^^^t'^^^ 
implies /igxt ^ 1, also (low- angle) domain theory is ap- 
pHcable up to the cross-over to Scenario II. (As m2 tends 
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towards one in Scenario II solely domain theory is appli- 
cable, the low- angle approximation has to be dropped - 
in particular for the wall energy.) 

Let us mention another observation supporting the 
conjecture that anisotropy effects are most prominent 
close to the field strength where the concertina vanishes: 
For Q ^ ^~^t^ the ground state for vanishing external 
field 

^ext — is no longer given by the uniform magneti- 
zation m = (±1,0,0), but a Landau or Concertina-type 
pattern, see Figure [38l has lower energy. The optimal 
periods w of the two latter pattern are determined by a 
balance of the wall energy and the anisotropy energy in 
the closure domains, and scale as ~ Q~^l'^{ttYl^ up 
to a logarithm. Hence we expect that in this regime, the 
concertina does not switch to m = (—1, 0, 0), but evolves 
to the pattern in Figure [38l 




FIG. 38: Continuous transition from the concertina pattern 
via the Landau state to the reversed concertina. Note that 
the total length of the walls and the Zeeman energy do not 
change while the anisotropy energy is smaller in case of the 
Landau state. 

In fact, that type of transition of the concertina pat- 
tern can be observed in CoFeB samples - that posses a 
stronger (transversal) uniaxial anisotropy, cf. Figure [4Ql 



VII. CONCLUSION 

In this work, we addressed the concertina pattern in 
very elongated thin-film elements. We provided an expla- 
nation of the formation and the coarsening of this pattern 
as the external field is reduced from saturation. 

We identified a parameter regime in which the uni- 
form magnetization becomes unstable to an oscillatory 
buckling mode. In this parameter regime, we derived a 
two-dimensional and thus numerically tractable reduced 
energy functional from three-dimensional micromagnet- 
ics. On the basis of the reduced model, we performed nu- 
merical bifurcation analysis: The bifurcation is slightly 
subcritical, but has a turning point, after which the buck- 
ling mode grows into the concertina pattern with its low- 
angle Neel walls. This is an alternative explanation for 



the formation of the concertina to the one proposed by 
van den Berg: An outgrow of an unstable mode instead 
of an in-grow of closure domains. Over a wide range 
of sample sizes, there is a good agreement between the 
explicit period of the unstable mode and the measured 
average period of the concertina pattern. In particular, 
the predicted dependence on film thickness and width is 
confirmed. However, the measured period exceeds the 
theoretically predicted one by a factor up to approxi- 
mately two. 

We gave an argument for this initial deviation that 
at the same time explains the coarsening: Domain the- 
ory based on the reduced model - where low-angle Neel 
walls are replaced by sharp discontinuity lines - shows 
that coarsened configurations are energetically favorable. 
More importantly, uncoarsened configurations eventually 
become unstable because the energy per period becomes 
concave. Based on the reduced model, we argued by a 
Bloch-wave Ansatz that this concavity indeed translates 
into a secondary instability of the concertina pattern 
with respect to long wave-length modulations. These sec- 
ondary instabilities are confirmed by numerical bifurca- 
tion analysis. The long wave-length instabilities are fur- 
ther confirmed by an extended bifurcation analysis that 
capitalizes on the near-degeneracy of the primary bifur- 
cation. This extended bifurcation analysis also showed 
that that the long wave-length instability of the primary 
branch extends all the way down to the turning point. 
Hence at the moment of its appearance, the concertina 
pattern already has a resulting period larger than the one 
of the unstable mode. That qualitatively explains the de- 
viation between the period of the unstable mode and the 
measured period of the concertina. Incidentally, these 
secondary instabilities are an asymmetric (with respect 
to the wave number) version of the Eckhaus instability 
introduced in the context of convective problems. 

We gave yet another argument for the deviation of 
the period of the unstable mode from the measured pe- 
riod of the concertina at its formation. Based on the 
reduced model, we established a continuous transition 
from the magnetization ripple, which is triggered by the 
polycrystalline structure of the material, and the con- 
certina pattern. On the level of the reduced model, the 
effect of an easy axis that varies from grain to grain trans- 
lates into a random transversal external field that smears 
out the subcritical bifurcation. Hence for a sufficiently 
strong ripple effect, as the concertina pattern becomes 
discernible from the ripple, it has already coarsened. 

Finally, we investigated the effects of a weak uniax- 
ial material anisotropy on the concertina pattern. We 
distinguished three effects: 1) a shift of the critical field 
that changes its sign already for weak anisotropics, 2) a 
change in the coarsened concertina pattern from "limited 
by shape anisotropy" to "limited by material anisotropy" 
that kicks in for somewhat larger anisotropics, 3) a 
change from a subcritical to a supercritical bifurcation 
for a sufficiently large transversal anisotropy. 
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A Scenario I 




FIG. 39: Table of scaling behavior of the optimal and marginally stable period and the amplitude of the transversal component 
in the regime ti'^ < Q < d'^^^r'^^H^^^ . 

The various analyses render a fairly complete picture 
of the energy landscape that in particular explains the 
hysteresis of the concertina pattern. 
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FIG. 40: Experiment: Hysteresis of a CoFeB sample of 60nm thickness and 30^im width. Following the coarsening we observe 
a transition to a Landau state at external field which turns into a concertina which degenerates and refines, and finally 
disappears. 
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